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Abstract. Quasi-periodic intensity disturbances propagating upward along the coro- 
nal structure have been extensively studied using EUV imaging observations from 
SOHO/EIT and TRACE. They were interpreted as either slow mode magnetoacous- 
4^ ■ tic waves or intermittent upflows. In this study we aim at demonstrating that time series 

of spectroscopic observations are critical to solve this puzzle. Propagating intensity and 
Doppler shift disturbances in fanlike coronal loops are analyzed in multiple wavelengths 
using sit- and- stare observations from Hinode/ElS. We find that the disturbances did not 
cause the blue- wing asymmetry of spectral profiles in the warm (~1.5 MK) coronal 
lines. The estimated small line-of- sight velocities also did not support the intermittent 
upflow interpretation. In the hot (~2 MK) coronal lines the disturbances did cause the 
blue- wing asymmetry, but the double fits revealed that a high- velocity minor component 
is steady and persistent, while the propagating intensity and Doppler shift disturbances 
are mainly due to variations of the core component, therefore, supporting the slow wave 
interpretation. However, the cause for blueward line asymmetries remains unclear. 



1. Introduction 

Upwardly propagating quasi-periodic i ntensity disturbances ar e known to exist in coro- 
nal loops (Bergh mans & Clettelll999l : I5e Moortel et al.ll2000h and have been studied 



extensively using EUV images from SOHO/EIT and TRACE (see a recent review by 



de Moorte l 2009). These propagating d isturbances were interpreted as slow magnetoa 



coustic waves dNakariakov et alJl2000u since their propagating speeds (~100 km s l ) 



are close to the sound speed in the corona and periodicities appear depending on the 
underlying driver such as the leakage of the photospheric p-modes (iDe Moortel et al 
20021 : IDe Pontieu et alJl2005h . 



However, similar phenomena were re-discovered with Hinode/XRT recently and 
interpreted as high-speed outflows because these distur bances were found associated 
with outflows at the loop's footpoints with Hinode/ElS (ISakao et al.ll2007l : lHarra et aT 



l2008h . Based on EIS observations of blueward asymmetries of spectral profiles, some 
authors suggested that these propagating dist urbances are most likely transient outflows 



produced by episodic coronal heatin g events dMcIntosh & De Pontieull2009l : lDe Pontieu et al 



2009; IDe Pontieu & Mclntoshl 12010) . How ever, the blueward sp ectral asymmetry is 



also consistent with the wave interpretation (iVerwichte et al]l201Qu . 
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Figure 1. (a) EIS Fe xn 195 A raster image, (b) Doppler velocity map. (c) Skew- 
ness map with a constant background subtracted. The contours shown in (a)-(c) 
indicate an area with the large negative skewness, and the thick vertical line marks 
the slit position for the sit-and-stare observations, (d) The normalized line profiles at 
positions A (in black) and B (in red) marked in (c). 



Wang et al.l f2009) first studied these propagating disturbances using the EIS time 



series and found that the Doppler shift variations have small amplitudes of only about 
1-2 km s" 1 , and are in phase with the intensity variations. These features are consistent 
with slow magnetoacoustic waves propagating upwards along coronal structures. In 
this study we examine the effect of the disturbances on the evolution of spectral profiles. 
The purpose is to determine whether the apparent propagating disturbances are the slow 
waves or high- velocity upflows, or possible coexistence of both. 



2. Observations and Data Analysis 

The observations were obtained on February 1, 2007 in AR 10940. The analysis of EIS 
sit -and-stare observat ions for propagating disturbances in Fe xn 195. 12 A was presented 
by lWang et al.1 (l2009h . They found quasi-periodic (12 and 25 mi n) propagating distur- 



bance s in both intensity and Doppler shift along the fanlike loops. iDe Pontieu & Mclntoshl 
(|2010h performed a "Red Minus Blue" (R-B) profile asymmetry analysis on the Fexn 
195. 12, Fe xm 202.04, and Fe xiv 274.20 A spectral lines, and discussed the cause of the 
quasi-periodicities in the R-B asymmetry using the forward modeling and R-B guided 
double fits techniques. Using a different approach, our study focuses on direct exam- 
ination of the influence of individual disturbances on the line profile shape. We also 
examine the line asymmetry using the skewness. The skewness is a typical measure 
of line asymmetry, defined as j [(A - Ao)/cr] 3 I(A)dA/ j I(A)dA, where Aq and cr are the 
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Figure 2. Propagating disturbances in coronal loops observed with EIS in the 
Fexn 195 A line, (a) Time series of relative intensity along the slit with a scale 
ranging from -10% (dark) to +10% (bright), (b) Time series of relative Doppler ve- 
locity with a scale ranging from -3 to +3 km s" 1 . (c) Time series of skewness with 
removal of a background which is taken as an average over time at each slit position. 
The scale range is from -0.3 to +0.3. (d) and (e) Comparisons of the normalized 
line profiles at positions PI and P2. The red line represents the spectral line profile 
at the bright disturbance, while the black line shows the profile in between (referred 
to as the background). 



mean and standard deviation of the line. FigureQ]shows the Fe xn raster image, Doppler 
shift, and skew ness maps. Since the Fexn 195.12 A line is self-blended with a weak 
line (195.18 Ai lYoung et al.ll2009h . a constant background has been removed from the 
skewness map. The skewness with large negative values indicates line profiles having 
evident blueward asymmetries (Figures [it and d). To find the change of a spectral pro- 
file in shape relative to the background profile, we first interpolate the spectral profile, 
and then normalize it to the peak intensity. We calibrate Doppler velocities measured 
with EIS in different emission lines by choosing Six 258.37 A in the long wavelength 
band and Fexn 195.12 A in the short wavelength band as th e reference lin e s and adopt- 
ing the relative wavelengths measured above the limb by 1 Warren etail (1201 lh . We 
adjust the wavelength scale of the reference lines by letting the quiet corona have no 
net velocities in these lines. 



3. Results for Spectral Profile Examination in Multiple Wavelengths 
3.1. Warm Coronal Lines: Fe xn 195 A and Fe xm 202 A 

Figure Q] shows that the strong blueward asymmetry of Fexn 195 A mainly occurred 
at the footpoints of weak loops associated with upflows, while in bright fanlike loops 
associated with downflows, signatures of the blueward asymmetry are less evident. Fig- 
ure [2] shows the propagating disturbances in intensity and Doppler shift of Fexn along 
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Figure 3. Comparison of the line profiles of Fexn 195.12 A on and off a bright 
propagating disturbance (PD) at position PI (top) and P2 (bottom) along the slit 
shown in Figure^, (a) Single Gaussian fit to the line profiles (red for PD, black 
for the background), (b) Same as (a), but for double Gaussian fit with the second 
component for a self-blended line, Fexn 195.18 A. (c) Triple Gaussian fit to the 
PD line profile (red line), by taking the background profile (fitted with the double 
Gaussian in solid black line) as a fixed component, where the excess component (in 
solid green line) is expected from the hypothetical upflow. (d)-(f) Same as (a)-(c), 
but for line profiles at position P2. In each panel, the reduced chi-squared statistics 
of the fits are marked. 



the bright fanlike loops. No propagating features found in time series of the skewness 
indicate that the disturbances did not cause changes in the line shape. This is confirmed 
by comparisons of the normalized line profiles at times on and off the propagating 
brightenings (Figures EJl and [2^). 

The line asymmetry also can be examined by the multiple fit analysis. Figure [3] 
shows that single Gaussian fits did not fully capture the data (with large reduced chi- 
squared statistics) because of the blending. The double Gaussian fits including the 
blended line of Fexn 195.18 A match well to the data. For such line profiles, it is 
impossible to add an extra component in fr ee multiple Gaussian fits. To check the sug- 
gestion by iDe Pontieu & Mclntoshl (1201 Oh that bright disturbances are caused by the 
quasi-periodic transient upflows, we make the hypothesis that the excess of the line 
profile for a bright disturbance (also called transient jets by some authors) relative to 
that taken before or after this disturbance at the same position (referred to as the back- 
ground) is fully contributed by an upflow, and then we can make a multiple Gaussian fit 
to the line profile for the disturbance by fixing the background components (Figures [3]: 
and[3]F). From the fits, we can estimate the line-of-sight velocity of the hypothetical 
upflow, which is found to be about 7 and 20 km s" 1 relative to the background com- 
ponent for the positions at the footpoint and the leg of the loops, respectively. We also 
notice that the Gaussian width of the fitted upflow component is larger than that of the 
background component by 16% at the position of the loop's footpoint, whereas it is 
smaller than that of the background component by 8% at the position of the leg. 

Although the Fexn 195 A observations have the highest signal-to-noise ratio, the 
self-blending made its analysis complicated as shown above. To confirm the results 
found in Fexn 195 A, we compare with the analysis of a clean line, Fexm 202.05 A. 
Figures and St* show that the propagating intensity and Doppler shift disturbances 
observed in Fe xm are very close to those seen in Fe xn. The negative skewness in the 
region y = 60 // -120 // indicates that the blueward line asymmetry may result from the 
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Figure 4. Time series of (a) relative intensity, (b) relative Doppler shift with the 
background trend subtracted, (c) Doppler shift, and (d) skewness in Fexm 202.05 A 
for the propagating disturbances in coronal loops observed with EIS. (e) Line profiles 
along with single Gaussian fits at PI (red line for PD, black line for the background). 
(f) Comparison between the normalized line profiles on and off the PD. (g) The 
PD line profile along with Double Gaussian fit (red line), by taking the background 
line profile (fitted with single Gaussian in black line) as a fixed component, where 
the excess component (green line) is expected from an hypothetical upflow. (h)-(j) 
Same as (e)-(g), but for position P2. 



associated upflows in this region (Figures St andHJl). However, the disturbances can 
propagate high above the region of the blueward asymmetry, and the skewness evolution 
shows no evidence for similar propagating features. These facts suggest that the upflows 
near the footpoints of the loop are persistent and stable during the observations, which 
cannot account for the quasi-periodic disturbances propagating to a high level. The 
line profile analysis shows that the estimated line-of-sight velocity of the hypothetical 
upflow is about 13 km s" 1 at the footpoint (Fig. |4^), while about 22 km s" 1 at the leg 
of the loop (Figure |4j). This result agrees well with that for Fexn. The line width of 
the flow component estimated for Fe xm is 42% larger at the footpoint, but 40% smaller 
at the leg of the loop than that of the background, which follows the same rule as for 
Fexn. In addition, the line-of-sight velocity of the hypothetical upflows estimated for 
Six 258. 37A (1.3 MK) are 11-13 km s" 1 , also consistent with the above results. 

Under the hypothesis that the propagating disturbances are transient upflows, the 
line-of-sight component of 10-20 km s" 1 along with the transverse component (as- 
sumed to be the observed propagating speed) of 100 km s" 1 implies an inclination of 
the magnetic field in c oronal loops to the vertical in the range 79°-84°, which is too 
large to be acceptable. iMarsh et all (l2009h measured the 3D geometry of fanlike loops 
using STEREO EUVI observations, and found their inclination to be 37° + 6° to the 
local normal. Moreover, the fact that the line width of the hypothetical upflow at the 
footpoint of the loop is much larger than that at the leg is unexpect ed, while this ma y 
be explained by the presence of multiple harmonics of the waves d Wan get al .1120091 ). 



232 



Wang, Ofman, and Davila 



(a) Fe XV 284.1597 % (b) Relative Doppler velocity km/s (c) Doppler velocity km/s (d) 
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Figure 5 . Same as Fig. 01 but for the Fe xv 284. 1 6 A line. 



Therefore, the above upflow hypothesis is not supported, instead the wave interpretation 
is preferred. 

3.2. Hot Coronal lines: Fe xiv 274 A and Fe xv 284 A 

The propagating disturbances along hot coronal loops in the Fe xiv 274.20 A and Fe xv 
284.16 A emission lines look similar, but appear different from those seen in the warm 
lines discussed in the last section. We find that the propagating disturbances detected 
in intensity and Doppler shift between warm and hot coronal lines are less correlated 
(not shown). Figures [5^ and 03 show that the propagating disturbances are weak in 
intensity, but more evident in Doppler shift and located near the footpoints of the loop. 
The presence of persistent upflows is revealed in these hot coronal lines, which are 
characterized by evident blueward line asymmetries in Fe xv indicated by the skewness 
(Figures[5]: and0i). Single Gaussian fits show evident excess emission in the blue wing 
above 100 km s" 1 (Figures [5^ andOi). Comparisons of the normalized line profiles 
reveal the enhancement in the blue wing of Fe xv due to the disturbances (Figures [5]F 
and[3). But this feature was not found in the warm coronal lines (Figures |4j and 0). 

Using the double fits by taking the background line profile as a fixed component, 
we estimate the line-of-sight velocity of the hypothetical upflow to be 97 and 78 km s" 1 
for two selected disturbances (Figures 0* and [5)). However, it is noticed that the line 
profiles for both the disturbance and the background have evident blueward line asym- 
metries, indicating the prese nce of a persi stent upflow. We use the constrained double 
Gaussian fits as suggested in IPeterl (l2010h to analyze the effect of the disturbances on 



the evolution of the core component and the second minor component. Figure [6| shows 
the line profiles of Fe xv 284 A along with the double fits upon and after the distur- 
bances, where the second minor component was forced to have a blueshift of at least 
70 km s" 1 . We find that the intensity increase of the core component accounts for about 
75% of the total increase in Fexv intensity for case PI, while it is more than the to- 
tal increase for case P2. The core component was blueshifted by 4 km s" 1 due to the 
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(a) Upon a propagating disturbance (b) After a propagating disturbance 
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Figure 6. Spectral profiles of Fe xv 284.16 A recorded in one spatial pixel at po- 
sition PI along the slit (a) upon and (b) after a propagating disturbance as shown in 
FigureEl (c)-(d) Same as (a)-(b), but for position P2. The constrained double Gaus- 
sian fits were made by forcing the minor component (red line) to have a blueshift of 
at least 70 km s" 1 . Vertical dashed lines indicate the peak positions of the two fitted 
components. The measured intensity (I), line width (W), and Doppler shift (V) for 
the fits are marked in each panel. 

bright disturbance in both cases, while the Doppler velocity of the minor component 
kept at about 100-120 km s" 1 with a small decrease upon the disturbance relative to 
the background (taken after the disturbances). This result indicates that the propagating 
intensity and Doppler shift disturbances with small amplitudes obtained in single fits 
are mainly ascribed to variations of the core component (loop emission) but not of the 
second component (high-speed upflows) in double fits for the hot coronal lines. There- 
fore, our result supports the interpretation of the propagating disturbances as slow mode 
magnetoacoustic waves propagating in coronal loops. 



4. Discussion and Conclusion 

We have presented a detailed analysis of propagating intensity and Doppler shift distur- 
bances in fanlike coronal loops using EIS sit-and-stare multi-wavelength observations. 
We adopted several techniques such as the skewness, the line profile normalization, and 
the multiple Gaussian fits to examine changes of the spectral profile shape caused by 
these disturbances in warm coronal lines (e.g., Fexn 195 A and Fexm 202 A) and in 
hot coronal lines (e.g., Fexv 284 A). 

Blue-wing asymmetries associated with Doppler blueshifts are found to exist near 
the footpoints of coronal loops from the skewness analysis, they are more evident in 
hot lines. No propagating features found in time series of the skewness suggest that 
the associated upflows are persistent during the observations. We find that these distur- 
bances appear uncorrected between the warm and hot lines, suggesting that they may 
propagate in different loops of different temperatures. The warm coronal lines show no 
evidence that the disturbances led to blue-wing asymmetries. The observed high (100- 
120 km s" 1 ) propagating speeds are not likely to be the result of hypothetical periodic 
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upflows whose line-of-sight velocity is estimated to be only about 10-20 km s" 1 . The 
hot coronal lines have a high- velocity (~100 km s" 1 ) component persistently existing 
in the blue wing during the period of observations, suggesting the presence of a steady 
upflow in hot coronal loops. The restricted double fits show that the propagating inten- 
sity enhancements are mainly contributed from the increase in the core component but 
not due to variations of the second minor component. Variations of the core component 
in intensity and Doppler shift having small amplitudes are consistent with propagating 
slow-mode waves. However, the disturbance-related excess in the blue wing of hot 
coronal lines revealed by comparisons of the normalized line profiles (Figures [3* and 
0) could not be explained by the superposition of the core component with Doppler 
shift oscillations caused by a slow wave and the background component produced by 
a steady upflow or static plasma based on a simple forward modeling, and the cause 
requires further investigation. 

Long-period (12 and 24 min) slow waves are not likely to be the direct result 
of photospheric p-mode leakage. Instead their association with steady upflows may 
suggest that these waves are produced by the same mechanism that drives the steady 
upflo ws. As these upf lows are progressively larger in lines formed at higher tempera- 
tures, [DelZanna (2008) suggested that they could be due to chromospheric evaporation 



following magnetic reconnection. The generatio n of slow waves is expected since they 
are a natural response to the impulsive heating (iNakariakov et al.1 [20041 : iTaroyan et al 



l2007h . Therefore, observations of these waves may pr ovide a potential diagnostic of 



coronal heating. More recently, iDel Zanna et al.1 (1201 lh proposed a rarefaction wave 



model based on reconnection between field lines of active regions and the nearby coro- 
nal hole, which can also explain the existence of persistent upflows and slow waves. 
But the reason why the upflows are persistent while the waves appear as quasi-periodic 
pulses is unknown and should be further studied in the future. 
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